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ABSTRACT 



^ : 

Context. On the nature, redshift, stellar populations and dust properties of optically faint or non-detected extremely red objects. 
Q ■ Aims. We determine the nature, redshift, stellar populations and dust properties of optically faint or non-detected extremely red objects (ERO) 
S-H ' found from our survey of the lensing clusters A1835 and AC114 (Richard et al. 2006). Comparison with properties of related galaxies, such as 
C/3 , IRAC selected EROs and a z ~ 6.5 post-starburst galaxy candidate from the Hubble Ultra Deep Field. 

• Methods. Using an updated version of Hyperz (Bolzonella et al. 2000) and a large number of spectral templates we perform broad-band SED 
^ fitting. The photometric observations, taken from Hempel et al. (2007), include deep optical, ACS/HST, ISAAC/VLT, IRAC/Spitzer data, and 
• »"H . for some objects also 24 pm MlPS/Spitzer and sub-mm data. 

' Results. For most of the lensed EROs we find photometric redshifts showing a strong degeneracy between "low-z" (z ~ 1-3) and high-z (z ~ 
| 6-7). Although formally best fits are often found at high-z, their resulting bright absolute magnitudes, the number density of these objects, and 
in some cases Spitzer photometry or longer wavelength observations, suggest strongly that all of these objects are at "low-z". The majority 
of these objects are best fitted with relatively young (S 0.5-0.7 Gyr) and dusty starbursts. Three of our objects show indications for strong 
extinction, with A v ~ 2A-4. The typical stellar masses of our objects are M t ~ (0.5 - 5) x 10 10 M Q after correction for lensing; for the most 
extreme ERO in our sample, the sub-mm galaxy SMMJ14009+0252 most likely at zm ~ 3, we estimate ~ 6. x 10 U M Q . For dusty objects 
star formation rates (SFR) have been estimated from the bolometric luminosity determined after fitting of semi-empirical starburst, ERO, and 
ULIRG templates. Typically we find SFR ~ (1 - 18) M Q yr 1 . Again, SMMJ14009+0252 stands out as a LIRG with SFR ~ 1000 M Q yr 1 . 
Finally, we predict the mid-IR to sub-mm SED of the dusty objects for comparison with future observations with APEX, Herschel, and ALMA. 
Concerning the comparison objects, we argue that the massive post-starburst z ~ 6.5 galaxy candidate HUDF-J2 showing observed properties 
very similar to our EROs, is more likely a dusty starburst at z ~ 2.3-2.6. This interpretation also naturally explains the observed 24 pm emission 
from this object and we predict its IR to sub-mm SED. 

Both empirically and from our SED fits we find that the IRAC selectee EROs from Yan et al. (2004) show very similar properties to our lensed 
EROs. Reasonable fits are found for most of them with relatively young and dusty stellar populations. 
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1. Introduction 

Various searches for distant galaxies based on a combination of 
deep broad-band optical and near-IR imaging using the Lyman 
break technique have been undertaken during the last years, 
e.g. based on data in the Hubble Ultra Deep Field (HUDF), on 
the GOODS survey, or others (e.g. Stanway et al. 2003, Yan et 
al. 2003, Bouwens et al. 2006). Relying on the use of strong 
gravitational lensing provided by rich foreground galaxy clus- 
ters, we have recently undertaken such a pilot program with 
the main aim of identifying z ~ 6-10 star forming galaxies 
(see e.g. Pello et al. 2004, Richard et al. 2006, and an overview 
in Schaerer et al. 2006). Such candidates are selected through 
the now classical Lyman break technique as optical drop-out 
galaxies showing an intrinsically blue UV restframe colour, as 
measured from near-IR colours. 

As a "by-product" galaxies with red spectral energy distri- 
butions (SEDs) are also found among the drop-outs. For ex- 
ample, in the study of two lensing clusters A1835 and AO 14 
by Richard et al. (2006) we found 8 lensed galaxies with 
R-Ks > 5.6 and red near-IR colours satisfying one of the often 
used criteria of "Extremely Red Objects" or EROs. The present 
paper focuses on these galaxies and related objects, using new 
ACS/HST observations in the zssolp band and Spitzer imag- 
ing obtained recently and discussed in Hempel et al. (2007, 
hereafter paper I). Down to the available depth (/ab ~ 27.3 to 
Vab ~ 28.) none except one of these objects is detected short- 
ward of zg50LP (Ave ~ 9100 A). This could imply that some 
of them are red galaxies at very high redshift (z <: 6) or lower 
redshift objects with a very strong extinction. Quantifying the 
properties of these EROs and comparing them with similar ob- 
jects is the main aim of the present work. 

Generally speaking, at z ^ 1 EROs are found be to either 
dusty starbursts or old passive galaxies. They are interesting in 
their own right and in the context of galaxy formation and evo- 
lution (e.g. review by McCarthy 2004). While normally EROs 
are found at relatively low redshift (z ~ 1-2), there are at- 
tempts to search for similar galaxies at higher z or for even 
more extreme - i.e. redder - objects, by selection at longer 
wavelengths. E.g. Yan et al. (2004) have identified IRAC se- 
lected EROs (or IEROs) in the HUDF. In these ultra-deep im- 
ages the IEROs turn out to be very faint in optical bands (~ 
27-30 mag in V, i or zssolp), such that they could be taken for 
optical drop-outs in images with less depth. Their redshift has 
been estimated between ~ 1.6 and 2.9 (Yan et al.). Other selec- 
tion criteria, such as searches for optical drop-out objects with 
old populations yield at least partial overlap with EROs, as e.g. 
demonstrated by the post-starburst z ~ 6.5 galaxy candidate of 
Mobasher et al. (2005), which also fares among the IEROs just 
mentioned. Also, a fraction of the sub-mm galaxies, detected 
through their strong dust emission, show optical to near-IR 
colours compatible with EROs (cf. e.g. Blain et al. 2002). This 
is also the case for one of our objects, SMMJ14009+0252, a 
known lensed sub-mm galaxy which also classifies as an ERO. 
Finally, optically faint and red objects have also been found by 
selection of high X-ray to optical fluxes. These so-called EXOs 
are thought to be AGN at very high redshift (z ^ 6) or in very 



dusty and/or sub-luminous host galaxies at more moderate red- 
shifts (z ~ 2 - 3; see e.g. Koekemoeret al. 2004). 

Given this variety of optically faint or undetected and red 
galaxies it is of interest to determine and compare their proper- 
ties to clarify their nature, and ultimately to obtain a coherent 
picture of these seemingly different galaxy populations. With 
these objectives, we have carried out a detailed quantitative 
analysis of the stellar populations and dust properties of the 
lensed EROs found by Richard et al. (2006) benefiting from 
the new ACS/HST and Spitzer photometry available for these 
lensing clusters (see paper I). The same SED fitting method, 
based on the photometric redshift code Hyperz (Bolzonella et 
al. 2000) and using a large set of spectral templates, was also 
applied to objects with similar SEDs, such as the IRAC selected 
EROs of Yan et al. (2004) and other objects from the HUDF. 

The paper is structured as follows. In Sect. [2] we briefly 
summarise the observational data for our lensed EROs. The 
SED fitting method is described in Sect. [3] Results for the 
SMMJ14009+0252 galaxy are presented and discussed in Sect. 
|U The other EROs are discussed in Sect. [5] In Sect. [6] we 
analyse and discuss the properties of related objects from the 
HUDF. Our main conclusions are summarised in Sect. [8] 

Throughout this paper we adopted the following cosmol- 
ogy: Q m = 0.3, H o =70 km s^'Mpc -1 in a flat universe. Except 
if mentioned otherwise, all magnitudes are given in the Vega 
system. 

2. Observations 

The present work deals primarily with the 8 EROs detected 
as red (R - Ks > 5.6) optically non-detected objects in the 
fields of the lensing clusters A1835 and AO 14 by Richard et 
al. (2006). The images obtained by Richard et al. were com- 
plemented with new ACS/HST observations in the zssolp band 
(optical), and where possible with Spitzer observations at near- 
to mid-IR wavelengths (IRAC 3.6, 4.5, 5.8, and 8.0 jum, and 
MIPS 24 fim photometry). 

All details concerning the photometric data reduction are 
given in Hempel et al. (2007, paper I). The photometry from 
paper I is summarised in Tables Q] [2] and [3] In this paper 
the EROs from Richard et al. were reselected from the Ks 
band image and optical and near-IR photometry measured with 
SExtractor using AUT0_MAG. For the photometry from differ- 
ent instruments source matching was done using object coordi- 
nates based on astrometry using the ESO-USNO-A2.0 catalog. 
For Spitzer total aperture-corrected flux densities were deter- 
mined from measurements in 3 " diameter apertures. No at- 
tempt has been made to correct the differences in flux measure- 
ment that are caused by the different apertures and image qual- 
ity. This is not necessary as we do not use aperture photometry 
for the optical and near-infrared images, or in case of Spitzer 
perform aperture correction. The photometry used, AUTO _MAG, 
is determined by measuring the flux in a flexible elliptical aper- 
ture around each object and accounts for the extended bright- 
eness distribution for the brighter objects. Possible differences 
in the photometry with respect to the earlier measurements of 
Richard et al. (2006) are discussed in paper I. Using fixed aper- 
ture near-IR photometry implies only relatively small changes 
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Table 1. Photometry of the selected ERO subsample in Abell 1835 and AC 114 taken from Hempel et al. (2007, paper I). All 
magnitudes are given in the Vega system. For conversion to the AB system see the filter properties listed in Table [3] Non- 
detections (lower limits) are 1 cr values. NA stands for non available. 



Object 


V 


R 


I 


Z850LP 


SZ 


J 


H 


Ks 


A1835-#l 


>28.1 


>27.8 


>26.7 


25.70±0.07 


24.44±0.27 


22.76+0.16 


22.40+0.08 


20.74+0.02 


A1835-#2 


>28.1 


>27.8 


>26.7 


>27.46 


24.08±0.26 


24.41 + 1.34 


21.78+0.06 


20.45+0.02 


A1835-#3 


>28.1 


>27.8 


>26.7 


24.06±0.07" 


23.78±0.10 


24.32+0.47 


22.55+0.07 


21.58+0.03 


A1835-#4 


>28.1 


>27.8 


>26.7 


25.48±0.14 


24.44+0.15 


23.56+0.18 


22.90+0.07 


21.95+0.03 


A1835-#10 


>28.1 


>27.8 


>26.7 


25.56+0.11 


24.00+0.12 


23.72+0.26 


23.36+0.13 


21.67+0.03 


A1835-#ll 


>28.1 


>27.8 


>26.7 


>27.46 


>26.9 


23.92+0.37 


23.49+0.18 


21.29+0.03 


A1835-#17 


>28.1 


>27.8 


>26.7 


>27.46 


>26.9 


> 25.6 


23.51+0.16 


22.11+0.03 


AC114-#1 


>28.5 


>27.7 


>26.8 


24.55±0.07 


NA 


21.26+0.04 


19.75+0.01 


18.62+0.001 



" This source appears double in the zssolp image. The fainter component has 25.31+0.07. In the SED modeling the brighter magnitude or the 
sum of two have been used. 



Table 2. IRAC and MIPS 24 fim photometry of the EROs listed in TableQ~|taken from paper I. All fluxes are given in /iJy. Upper 
limits are 1-sigma noise values at the position of the sources. Due to source blending the data is uncomplete, i.e. not available for 
objects A1835-#3, A1835-#10, and A1835-#l 1, and partially for A1835-#l. 











IRAC/Spitzer 








MIPS 


Object 




3.6 


4.5 




5.8 




i 


;.o 


24.0 


A1835-#l 


2.9 = 


t 0.2 


4.0 


±0.2 




blended 




< 10 


A1835-#2 


14.4 = 


t0.2 


23.; 


5 ±0.3 


37.6 ± 


1.5 


50.9 


± 1.6 


320 ± 1 1 


A1835-#4 


1.9 = 


t0.2 


1.0 


±0.3 


<0.7 




< 


0.6 


< 10 


A1835-#17 


2.4 = 


t0.2 


1.6 


±0.2 


< 0.6 




< 


0.5 


< 10 


AC114-#1 


67.3 = 


h0.5 


64.4 ± 0.5 


50.3 ± 


1.5 


44.9 


±2.2 


189.0 ±8.9 



and does not alter our overall conclusions, as test computations 
have shown. The new measurements, adapted to the morphol- 
ogy of the EROs, are used in the present paper. In the SED 
modeling discussed below we consider also a minimum photo- 
metric error to account for uncertainties due to matching pho- 
tometry from different instruments. 

No object from the subsample of "optical dropout" EROs 
discussed in paper I is detected shortward of the R band. For 
modeling in the present work we include the V band non- 
detection as the dropout constraint. Non-detections at shorter 
wavelengths are redundant and are therefore not included in 
the SED fitting. 

An overview of the observed optical, near-IR, and 
IRAC/Spitzer fluxes of all the objects discussed in this paper 
is shown in Fig.Q] 



3. SED fitting method 

An SED fitting technique based on an updated version of the 
Hyperz code from Bolzonella et al. (2000) is used to constrain 
the redshift, stellar population properties (age, star formation 
history), and extinction of the galaxies studied in this paper. 
To do so we closely follow the procedures outlined in Schaerer 
& Pello (2005). In addition we have included other synthetic, 
empirical and semi-empirical spectral templates, as described 
below. 



3.1. Photometry 

Ground-based, HST, and Spitzer photometry of the lensed 
EROs is taken from paper I. The following bands have been 
included in the SED fitting: V, R, I nondetections, zssolp from 
ACS/HST, SZ, J, H, Ks from ISAAC/VLT, channels \-A from 
IRAC/Spitzer and 24 fim MlPS/Spitzer where available. The 
filter properties are listed in Table [3] 

For other objects, included here for comparison, the origi- 
nal photometry was taken from the literature. 

In some cases a prescribed minimum photometric error is 
assumed, to examine the influence of possibly underestimated 
error bars and to account for uncertainties in absolute flux cal- 
ibrations when combining photometry from different instru- 
ments. 

3.2. Spectral templates 

The following spectral templates, assembled into several 
groups, have been used in the fit procedure. 

- Bruzual & Chariot plus Coleman et al. (1980) empiri- 
cal templates galaxies of all Hubble types (hereafter named 
BC or BCCWW group). The theoretical Bruzual & Chariot 
(2001) models, taken here for solar metallicity, include 
various star formation histories representative of different 
Hubble types (burst, and e-folding times of t=1, 2, 3, 5, 
15, 30, and oo Gyr corresponding to E, SO, Sa, Sb, Sc, Sd, 
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Table 3. Properties of the photometric filters used for the SED 
fitting of objects in the field of A1835 and of AC114. For 
AO 14 alternate filters are listed in parenthesis and SZ obser- 
vations are not available (entry NA). Col. 1 indicates the fil- 
ter name, col. 2 the effective wavelength in micron, col. 3 the 
effective bandpass (filter "width") in micron computed with 
a Gaussian approximation. AB corrections (Cab), defined by 
rriAB = ffivega + Cab, are listed in col. 3. Detailed information on 
the photometry can be found in paper I. 



Filter 




A/1 [/im] 


Cab 


V 


0.543 


0.056 


0.018 


R (^702) 


0.664 (0.700) 


0.075 (0.123) 


0.246 (0.299) 


/ (/ 8 i4) 


0.817 (0.807) 


0.117 (0.137) 


0.462 (0.445) 


Z850LP 


0.911 


0.114 


0.540 


SZ (NA) 


1.070 


0.094 


0.698 


J 


1.259 


0.167 


0.945 


H 


1.656 


0.180 


1.412 


Ks 


2.165 


0.181 


1.871 


IRAC 3.6 yum 


3.577 


0.427 


2.790 


IRAC 4.5 yum 


4.530 


0.567 


3.249 


IRAC 5.8 yum 


5.788 


0.801 


3.737 


IRAC 8.0 yum 


8.045 


1.634 


4.392 



and Im types). The IMF adopted in these models is the 
Miller-Scalo IMF from 0.1 to 125 M Q . We have not in- 
cluded the templates from the Bruzual & Chariot (2003) 
update, as this concerns mostly high resolution spectral li- 
braries which have no impact on our results. 

- Starburst SEDs from Schaerer (2002, 2003) models at 
different metallicities extended up to ages of 1 Gyr and 
considering instantaneous bursts or constant star formation 
(s04gyr group). The overall SEDs predicted by these syn- 
thesis models, including in particular nebular continuum 
emission neglected in the Bruzual & Chariot models, are 
basically identical to the ones from Starburst99 (Leitherer 
et al. 1999). However, a larger variety of metallicities is in- 
cluded. These models assume a Salpeter IMF from 1 to 100 
Mo. 

- Maraston models (Maraston 2005) including a semi- 
empirical treatment of thermally pulsating AGB stars, 
whose contribution in the near-IR may be significant in 
certain phases. This non-standard approach may therefore 
lead to different age and stellar mass estimates than other 
synthesis codes (see Maraston 2005 details). The templates 
used here include simple stellar populations (bursts) with 
ages up to 15 Gyr, and exponentially decreasing star for- 
mation histories with e-folding times up to 2 Gyr. 

- Empirical or semi-empirical starburst, ULIRG and 
QSO templates. In addition to starburst templates from the 
Calzetti et al. (1994) and Kinney et al. (1996) atlas included 
in the public Hyperz version, we have added the HST QSO 
template of Zheng et al. (1997), and templates of metal- 
poor H 11 galaxies SBS0335-052 and Tol 1914-266 includ- 
ing numerous strong emission lines (Izotov private commu- 
nication). 

To include also more obscured objects we have added 
UV to millimeter band templates of EROs, ULIRGS, star- 



burst and normal galaxies (HR 10, Arp 220, M82, NGC 
6090, M51, Ml 00, NGC 6949) from fits of GRASIL mod- 
els to multi-wavelength observations (Silva et al. 1998, 
named GRASIL group). These templates are therefore 
semi-empirical templates. This template group will be used 
in particular to predict mid-IR to sub-mm fluxes, and hence 
to estimated total bolometric luminosities, after fitting the 
optical to 8 /im part of the spectrum. 

3.3. Fit procedure and determined parameters 

The main free parameters we consider are: the spectral tem- 
plate (among a group), redshift z , and (additional) extinction 
(Ay) assuming a Calzetti et al. (2000) law. In our standard cal- 
culations z is varied from to 10, and A v from to 4 mag. For 
test purposes higher values of A v and other extinction laws are 
also allowed. 

To increase the diversity of empirical or semi-empirical 
templates and to allow for possible deviations from them, red- 
dening is optionally also considered as a free parameter. In 
this case, this obviously corresponds to an additional redden- 
ing. Test computations have shown an excellent consistency 
between photometric and spectroscopic redshifts using e.g. the 
GRASIL template group to fit the near-IR to IRAC observa- 
tions of Stern et al. (2006) of the ERO HR10, once allowing 
for possible additional reddening. A similar approach with em- 
pirical templates was also adopted by Rigby et al. (2005). As 
for all templates the corresponding dust emission is not treated 
consistently. 

Finally, from the luminosity distance of the object the scal- 
ing of the template SED to the observed absolute fluxes yields 
an absolute scaling property, such as the stellar mass or the star 
formation rate (SFR) when templates generated by evolution- 
ary synthesis models are used. 

To estimate stellar masses we use two different approaches. 
1) We determine the restframe absolute K band magnitude 
M rest (/Ti) from Hyperz and assume a typical light-to-mass ratio 
(Lk/M). The stellar mass is then determined as: 

M+ = 10- a4[A ™ s) - 3 - 3] /(WM) (1) 

in solar units. The numerical value 3.3 is the solar absolute Ks 
band magnitude. This approach is applicable to any template, 
theoretical or empirical ones. For comparison with other dusty 
galaxies, we adopt the value Lk/M = 3.2 used for SCUBA 
galaxies by Borys et al. (2005). 

2) When using templates from evolutionary synthesis mod- 
els, the stellar mass (and/or SFR) can be determined from the 
absolute scaling of the best fit template to the observed fluxes 
and to the best fit redshift. Note that the Bruzual & Chariot 
models used here assume a Miller-Scalo IMF from 0.1 to 125 
M , whereas the S04 models assume a Salpeter IMF from 1 to 
100 M Q . 

The star formation rate is a natural quantity when star for- 
mation over a certain time scale (or at a constant rate) is consid- 
ered. This quantity can therefore only be determined through 
SED fits using theoretical templates assuming constant star for- 
mation. Alternatively, for objects with good fits using the multi- 
wavelength GRASIL templates covering the restframe UV to 
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sub-mm domain, we will determine SFRs from the total bolo- 
metric luminosity derived over the available spectral range and 
applying a standard Kennicutt (1998) relation between Lboi an d 
SFR. 

Absolute quantities such as the stellar mass, SFR, and bolo- 
metric luminosity depending on the luminosity distance must 
also be corrected for the effects of gravitational lensing. The 
magnification factor of each source was determined using the 
mass models of A1835 (similar to Smith et al. 2005) and 
AC114 (Natarajan et al. 1998, Campusano et al. 2001), fol- 
lowing the same procedure as in Richard et al. (2006). Because 
of the slight dependence of the magnification with the source 
redshift z s , at the location of the EROs, we computed different 
estimates assuming z, = 0.5, 1, 2, 3, as well as 7 for com- 
parison. The different values, given in Table|4]for each object, 
reflect the uncertainty in the magnification factor, the source 
redshift being the dominant source of error. In any case, for the 
bulk of these sources located somewhat away from the cluster 
center (see Figs. 10, 11 in Richard et al.), the magnification is 
relatively small. 

Table 4. Magnification factors jj from the lensing models of 
A1835 and AC114 predicted for various source redshifts z s - 
The values of p. are dimensionless magnification factors, and 
not in magnitudes. 



Object 


z s = 0.5 


z s = 1. 


z s = 2. 


z s = 3. 


z s = 7. 


A1835-#l 


1.13 


1.21 


1.26 


1.28 


1.33 


A1835-#2 


1.38 


1.71 


1.95 


2.05 


2.32 


A1835-#3 


1.24 


1.43 


1.54 


1.59 


1.72 


A1835-#4 


1.25 


1.45 


1.57 


1.62 


1.76 


A1835-#10 


1.19 


1.33 


1.42 


1.45 


1.54 


A1835-#ll 


1.14 


1.23 


1.28 


1.30 


1.35 


A1835-#17 


1.14 


1.23 


1.28 


1.30 


1.36 


AC114-#1 


1.49 


2.24 


2.87 


3.15 


4.03 



4. The sub-mm galaxy A1835 - #2 
(SMMJ1 4009+0252) 

As already mentioned in Richard et al. (2006), this ERO cor- 
responds to the known sub-mm source SMMJ 14009+0252 
(Ivison et al. 2000, Smail et al. 2002, Frayer et al. 2004). With 
AC1 14-#1 this object is the brightest optical dropout ERO from 
our sample. 

4.1. SED fitting results 

The observed SED shows, in F v units (cf. Fig. HI, a continu- 
ously increasing SED from the near-IR, through the 4 IRAC 
channels and up to 24 /mi (MIPS), where this object is also 
detected. When all IRAC bands are included in fits, the solu- 
tions with the "standard" templates are driven to high-z (z p hot 
between 6 and 8). However, a reasonably low x 2 is achieved 
for all redshifts z p hot ~ 3; see^f 2 map on Figj2] For all redshifts 
a very large extinction (Ay ~ 3-4-I) is found as a best fit. 

As shown in Figs. [3] and |H the high-z solution (zfit=7.46) 
provides an excellent fit to the observed SED. This template 



corresponds to a young burst (6 Myr) + high A v . The best fit 
with GRASIL templates is obtained at z fit =2.78 with the NGC 
6090 template plus additional extinction of A v =1.4. Imposing 
a maximum redshift of 4 to the BCCWW templates, one finds 
a very similar best photometric redshift (zfn=2.95) for an el- 
liptical with 0.36 Gyr plus 2.4 mag extinction in A v . These two 
z ~ 3 solutions are also plotted in Fig. [3] showing a discrepancy 
at ~ 0.95-1.1 fjxa (cf. below). 

Actually the overall SED of this object, including in par- 
ticular our MIPS 24 /mi and the SCUBA measurements from 
Ivison et al. (2000), is rather well fitted with semi-empirical 
templates from GRASIL for redshifts z ~ 2.8-3., as shown on 
Fig. H Templates with very strong dust emission such as Arp 
220 are needed to reproduce the observed ratio of the sub-mm 
to near/mid-IR flux. For example, templates of more moderate 
starbursts like M82 and NGC 6090 underpredict the sub-mm 
emission. The Hyperz best fit with the Arp 220 template re- 
quires an additional extinction of A v = 1.4 (for the Calzetti 
et al. law). The only difficulty with fits at z&t~ 3 is the excess 
emission observed in the S Z band, which is ~ 5cr above the ex- 
pected level at such redshift. A natural explanation could be the 
Mg n ,12798 emission line seen in type 1 AGNs (e.g. Gavignaud 
et al. 2006). 

Already from the observed monotonous flux increase ob- 
served across all 4 IRAC bands it is quite clear that this source 
cannot be at redshift much smaller than ~ 3. Otherwise the typ- 
ical flux depression, associated to the transition from the stel- 
lar peak at 1.6 /mi (restframe) to the raising dust emission at 
longer wavelengths (cf. John 1998, Sawicki 2002) should be 
seen. Excluding higher redshift solutions from the near-IR to 
8 pm data used here for the Hyperz SED fitting is difficult, 
and would require a more complex stellar population plus dust 
modeling. However, if at z ~ 6-8 as suggested from the formal 
best fits, the absolute magnitude of this object would be rather 
exceptional (M res '(Ks) 29.6 or M res, (V) ~ -26.6 with- 
out correcting for lensing), rendering this case very unlikely. 
Furthermore radio and sub-mm data (cf. below) as well as our 
"global" SED analysis favour z <k 6. For these various reasons 
we conclude that the most likely redshift of this object is z ~ 
3. This redshift is larger than the estimate based on the radio- 
submm spectral index a^ 5 ®, but in agreement with the one from 
submm colours (cf. Ivison et al. 2000). It is also larger than 
our previous estimate based on optical to near-IR photometry 
(Richard et al. 2006). 

Assuming z ~ 3 and a magnification factor p = 2 (cf. Table 
lU one obtains the following estimates: with a rest-frame abso- 
lute magnitude Af est (Ks) ~ -27.7) this object is ~ 3.4 to 3.6 
magnitudes brighter than M* at this redshift (cf. Kashikawa et 
al. 2003) or than M*(Ks) from 2MASS in the local Universe 
(cf. Kochanek et al. 2001). 

The mass, estimated from the best fitting Bruzual & Chariot 
and S04 templates (with an age of ~ 0.36 and 0. 14 Gyr respec- 
tively), is M+ ~ 1.2 x 1O 12 //(M Q, slightly more massive e.g. 
than the most massive SCUBA galaxy discussed by Borys et al. 
(2005) and typically an order of magnitude more massive than 

1 For the value L K /M = 3.2 adopted for SCUBA galaxies (cf. Borys 
et al. 2005) one obtains M* ~ 7.8 x 10" /pM e 
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the most massive z ~ 3 Lyman break galaxy observed with 
Spitzer by Rigopoulou et al. (2006). 

Integrating the global SED of Arp 220 adjusted to the ob- 
servations (see Fig. |4]i and assuming z = 2.78 one obtains a 
total luminosity of Lboi ~ 1.2 x 10 13 /// L , close to the limit 
between ultra-luminous and hyper-luminous infrared galaxies 
(ULIRG and HyLIRG). Using standard SFR conversion fac- 
tors (Kennicutt 1998) this corresponds to an estimated star 
SFR ~ 1050 M yr-\ adopting yu ~ 2 (Table©. 

4.2. Discussion 

The radio-submm spectral index ff| 5 4 ° = 0.60 + 0.03 indicates 
a likely redshift of 0.7 < z < 2.3 (Smail et al. 2000, Ivison et 
al. 2000). From the 450- to 850-/im flux a coarse estimate of 
z ^ 2.8 has been derived by Hughes et al. (1998). From the 
near-IR to submm SED and from the low value of or? 5 P, Ivison 
et al. (2000) argue that A1835-#2 is more likely at 3 < z < 5 
and that the radio flux contains some AGN contribution. A re- 
cent reanalysis of the radio and submm flux by Aretxaga et 
al. (2003) yields a higher redshift estimate of z ~ 4.1 + 0.8 
However, their best fit SED strongly overpredicts the observed 
near- to mid-IR flux. Our earlier photometric redshift estimate 
of z ~ 1.2-1.6 (Richard et al. 2006) is now superseded by the 
present analysis including in particular the longer wavelength 
IRAC/Spitzer observations leading to a higher z. The various 
redshift estimates, including the one presented here, can be rec- 
onciled if the radio flux contains a contribution from an AGN, 
as already pointed out by Ivison et al. (2000). For our best red- 
shift, z ~ 3, the AGN contribution does not need to be strong, 
as also discussed by these authors. As already noted above, the 
observed SZ band excess could also be an indication for an 
AGN. 

At the high luminosities of this object, in the ULIRG range, 
the AGN fraction is high (> 40-50 %, cf. Veilleux et al. 1999, 
Alexander et al. 2004) rendering the AGN hypothesis quite 
likely. However, is there other direct evidence for an AGN ? 
Ivison et al. (2000) have obtained an upper limit in soft X-rays 
(0.1-2.0 keV) from ROSAT archival HRI observations. In our 
recent Chandra observations of Abell 1835, described in paper 
I, this object remains undetected with flux limits of the order 
of < (2. - 3.) x 10~ 16 erg s" 1 cm" 2 in the 0.5-7.0 keV band, 
for photon power law index T between 1.0 and 2.0. The corre- 
sponding limit for 2.0-10.0 keV and r = 1.4 is < 2.5 x 10~ 16 
erg s _I ctrT 2 . A comparison with the X-ray and 24 /im fluxes of 
starbursts and AGN compiled by Alonso-Herrero et al. (2004) 
places this X-ray limit well below the typical range of hard X- 
ray selected AGN. A more detailed analysis will be needed to 
examine how much room these new constraints leave for a pu- 
tative AGN in this object. 

Compared to other SCUBA galaxies studied also in the 
rest-frame optical (cf. Smail et al. 2004) A1835-#2 fares among 
the faintest ones in K and among the "reddest ones" in opti- 
cal/IR flux. With m{Ks) ~ 20.5 it is close to the faintest objects 
of Smail et al., which have 18 55 Ks 55 20.9; among the 7 
confirmed sub-mm galaxies observed observed in the SCUBA 
Cluster Lens Survey of Frayer et al. (2004) it is the second 



faintest object in K surpassed only by SMMJ00266+1708 with 
m(K) = 22.36 + 0.16, and the second reddest in J - K. After 
lensing correction the magnitude of A1835-#2 is m(/Ti) ~ 21.2. 
Several other sub-mm galaxies are known with very faint flux 
levels at K > 21.0-21.9 (Smail et al. 2002, Dannerbauer et 
al. 2002). Its restframe V-band to IR luminosity ratio is very 
low, ~ 5. x 10~ 4 , placing it among the 5 most extreme sub-mm 
galaxies when compared to the Smail et al. (2004) sample. 

In terms of stellar populations we find a dominant stellar 
age of ~ 0.36 Gyr or younger for A1835-#2, similar to the mean 
ages of ~ (3 10 - 530) + (80 - 90) Myr estimated by Smail et al. 
(2004) for a sample of sub-mm galaxies and optically faint ra- 
dio galaxies. The extinction we estimate (Ay ~ 2.4-3) is some- 
what larger than the average of Ay ~ (1.70-2.44)+(0. 13-0.14) 
found by Smail et al., but comparable to the median Ay ~ 
2.9 + 0.5 determined by Takata et al. (2006) for sub-mm galax- 
ies from the Balmer decrement. The best fit with the Arp 220 
spectrum, requiring an additional extinction of Ay = 1.4, also 
indicates that we're dealing with an object with a rather excep- 
tionally large extinction! 

We note also that the stellar mass estimated from the SED 
fit (M* ~ 1.2 X 1O 12 ///M ) is consistent with the mass being 
built up at the high SFR of ~ 2100/^ M yr~' over a period 
of 55 360 Myr. This leaves room for ~ 25 % of the stellar mass 
being formed from a previous star formation event. 

5. SED fitting results for other objects 

We now present the results from the SED fits for the individual 
objects. First we discuss in detail the objects for which Spitzer 
photometry (detections or upper limits) is available. The re- 
maining objects are addressed in Sect. 15.51 The main results 
are summarised in Table [5] and [6] The magnification factors /i 
needed to correct for gravitational lensing are listed in Table [4] 

5.1. A1835-#1 

Overall SED fits for this object are rather degenerate and of 
poor quality (high x 2 ), showing several minima for its photo- 
metric redshift, most of them at z 55 3 (see Fig. [5}. Formally, 
all template groups yield a x 1 minimum at znt ~ 0.4. However, 
the best fits are of lower quality (higher^ 2 ) than for the other 
objects, since the photometry yields an apparently somewhat 
"non-monotonic" SED. Furthermore the photometric redshift 
is only loosely constrained as no IRAC photometry is available 
for channels 3 and 4 (5.8 and 8.0 fim) due to blending with 
other sources. 

There is a strong degeneracy in this case between redshift 
and extinction (Pig. [5j. Two different solutions coexist, one at 
Zfit ~ 0.4 and Ay > 3.4, and another one at Zfa ~ 1.5 with 
Ay < 1.2. The first solution is strongly degenerate in the age- 
Ay plane as well, thus providing loose constraints on the stellar 
mass: the stellar ages vary quite strongly from ~ 2.3 Gyr (BC 
models) to ~ 10 Myr for the Maraston and S04gyr models, and 
the corresponding stellar masses range from ~ l.X lO IO /yuM 
(BC) to (0.4 - 1.1) x 10 8 lliM Q . Even if we consider the solu- 
tions around z ~ 1.5 there are significant uncertainties. E.g. the 
BC and Maraston models require little extinction (Ay ~ 0.6), 
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whereas S04gyr and GRASIL templates indicate a higher ex- 
tinction (even Ay ~ 3.8 for S04gyr!). Stellar ages of ~ 5.5 Gyr 
(1 Gyr) are found for the BC and Maraston (S04gyr) models; 
the corresponding stellar mass is estimated as ~ 8. x 1O 1O juM . 
For the S04gyr models one obtains M* ~ 6.3 x 10 11 fiM @ . 

For illustration we show several SED fits including with the 
semi-empirical GRASIL templates in Fig. [6] The latter allow us 
in particular to estimate the mid-IR to sub-mm flux. In particu- 
lar we note that the 24 fj.m non-detection probably rules out the 
very dusty solution at z ~ 0.4, rendering z ~ 1.5 more likely. 
However, for this object it is clear that the uncertainties on all 
derived parameters are large, and larger than for the other ob- 
jects discussed here. For this reason the entries in Table [5] are 
left blank for this object. For comparison we note that the em- 
pirical classification based on near-IR colours would indicate 
an "old passive" object (paper I). 

5.2. A1835-#4 

For this object the best fits are consistently found at low red- 
shifts, Zflt ~ 1.2 well constrained by the measurement of the 
stellar 1.6 fim peak measured in the IRAC channels. The cor- 
responding x 1 maps and the best fit SEDs for this object are 
shown in Figs.|7]and[8] 

Best fit templates correspond to bursts of ~ 4.5 Gyr with 
no extinction for Bruzual & Chariot models or to the ellipti- 
cal template from CWW, i.e. an old and dust-free galaxy. With 
the S04 templates the best fit is of similar quality, yielding a 
younger burst age (~ 0.6 Gyr) and some extinction (Ay ~ 1.6). 
From the BC and S04 model sets the estimated mass is M* ~ 
(1.1 - 1.7) x 1O 1O ///M , with the magnification factory ~ 1.5 
(cf. TableSJi. If we assume L K /M = 3.2 as for SCUBA galaxies 
(cf. Borys et al. 2005), one obtains M = 9.4 x 1O 9 ///M . For 
this object SED fits with Maraston models yield a solution of 
similar quality, but a lower redshift of z&t ~ 0.8. The other fit 
parameters a burst age of 1.7 Gyr, Ay = 1.2, and a stellar mass 
of 4.8 x lO 9 /yuM . 

For comparison the best fit to GRASIL templates is found 
at Zfl t =1.24 with an M82 template (and no additional extinc- 
tion), as shown in Fig. [8] Although the quality of this fit is less 
than the ones mentioned above, we cannot completely rule out 
the presence of dust. Observations at longer wavelengths would 
be needed. Assuming that the M82 template is valid, we esti- 
mate a bolometric luminosity of to Lboi ~ 2.6 x 10 10 /fiL G or a 
star formation rate of just S FR ~ 5 /fi M yr~ 1 . 

We note also that the results from the best fit agree with 
the empirical classification as "old passive" galaxy based on 
near-IR colours (see paper I). 

As a cautionary note, we remind the reader that this object 
has been found variable over ~ 1 month in the ISAAC photom- 
etry taken in the SZ band (see Richard et al. 2006). The SZ 
flux adopted here corresponds to the average between the 2 pe- 
riods. It is currently unclear if and to which extent the apparent 
variability influences the results derived here. 



5.3. A1835-#17 

This object is one of the few for which there is no ambiguity on 
the photometric redshift. See y 1 ma P on Fig|9] This is mostly 
due to the fact that the stellar 1.6 /urn peak is clearly observed 
between the Ks band and the first two IRAC channels (3.6 and 
4.5 /mi). With nearly all templates one obtains zm ~ 0.8; a 
somewhat lower value of zg t ~ 0.69 is obtained with templates 
from the S04gyr group. 

In all cases a very high extinction (Ay ~ 3^1) is needed. 
However, models with very different ages yield fits of similar 
quality (x 1 ): 7.5 Gyr with the BC models, 1 Gyr with Maraston 
models, and 10 Myr with the S04 templates. The SF histo- 
ries correspond to bursts in all of them. This age uncertainty 
is most likely due to the fact that we rely here on the pre- 
dictions in the rest-frame spectral range ~ 1.6-2.8 fim, where 
the evolutionary synthesis models are more uncertain than at 
shorter wavelengths. The stellar masses derived from the burst 
model fits are ~ (0.6 - 2.) x lO IO /yuM for the Maraston and 
BC templates and significantly smaller from the (younger) S04 
template (M* ~ 7. x 10 7 /^M Q ). For a value L K /M = 3.2 
adopted for SCUBA galaxies (cf. Borys et al. 2005) one ob- 
tains = 3.1 x 1O 9 M , much lower than the typical masses 
of SCUBA galaxies. The magnification factor for this object is 
li ~ 1.2 (cf. Table©. 

Using the semi-empirical GRASIL template group the best 
fits are found at z&t = 0.78 with the SED of the Sbc galaxy 
M 51 with an additional optical extinction of A v ~ 3.8. The 
corresponding mid-IR to sub-mm SED is shown in Fig. [TUl 
with a bolometric luminosity of 5.5 x 1O 9 //vL corresponding 
to SFR= 0.9//i M yr _1 . If this object is indeed a very dusty 
starburst, which can in principle be verified with longer wave- 
length observations, it is much fainter than SCUBA galaxies 
(SMG) at the same redshift — e.g. ~ 4 mag fainter in K (cf. 
Smail et al. 2004). 

5.4. AC114-#1 

Using the templates from synthesis models the best fits for this 
object are found between zm ~ 1 and 2.5, with a secondary 
though less likely solution at high redshift (see Fig. fTTT i. Over 
the interval ~ 1-2.5 the photometric redshift is actually not well 
determined, since the curvature of the SED measured in the 4 
IRAC bands is small and hence the position of the 1.6 fim peak 
- the main constraint on z - only loosely constrained. 

Both the BC and the Maraston templates give quite similar 
best fits: z& =1.3-1.5, a burst with a maximal age of ~ 3.5^4.5 
Gyr, a stellar mass of ~ (1.6 - 2.6) x 10 12 //jM q and a large 
extinction (Ay ~ 1.6-2.4). The magnification factor for this 
object is fi ~ 2.2. Best fits with a similar^ 2 are found using the 
S04gyr templates at Zfi t =1.6 for a burst of 0.9-1.0 Gyr age with 
Ay — 2.8 and a stellar mass corresponding to ~ 1 .3 x lO I2 //iM 
0. The main difference between these models is a lower age and 
higher extinction in the latter. 



2 A lower mass M t = 2.4 x x 10" ///Af is obtained using Eq. 1 and 
L K /M = 3.2. 
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Using the semi-empirical GRASIL templates yields a best 
fit at Zf\ t =0.9-1.0 for the M51 template. A strong additional 
extinction of Ay = 3.8 is required, and the overall fit is lower 
quality (higher x 2 ) than the fits discussed above. The overall 
SED resulting from these different fits is shown in Fig. Q~2] 
Interestingly the GRASIL template also reproduces quite well 
the observed MIPS 24 fim flux, although this was not included 
in the fit procedure. In any case the 24 /jm flux is a strong indi- 
cation for the presence of dust in this galaxy. In other words 
solutions with non negligible extinction at z ~ 0.9-1.5 are 
favoured by this additional constraint. 

If at z = 0.97 the GRASIL template shown in Fig. [HI has 
a bolometric luminosity of 2.8 x 10" //xLq (close to the LIRG 
range) corresponding to SFR= 48/^M , with /i ~ 2 (Table|4]l. 

From the GRASIL SED shown here we may also expect 
a fairly strong sub-mm flux, in the range detectable with cur- 
rent instrumentation. To the best of our knowledge the south- 
ern cluster AC1 14 has so far not been observed in this spectral 
range. 

5.5. Other EROs in A1835 - galaxies without Spitzer 
photometry 

The objects treated here are those from Table [6] for which con- 
tamination by neighbouring sources does not allow us to deter- 
mine photometry from the Spitzer images, i.e. the objects #3, 
#10, and #1 1 in Abell 1835. The main properties estimated for 
#3 and #10 are summarised in Table[6] 

The SED of #11, shown in Fig. [T] precludes any reliable 
photometric redshift estimate; above z ^ 1 good fits can be 
found at all redshifts. For this reason this object is not discussed 
further. 

As for other objects discussed earlier, #3 and #10 show a 
degeneracy between low and high-z, with x 1 minima found at 
Zflt ~ 1-1.5 and 5-6.0 for #10 (5.-6.5 for #3). However, if at 
high-z their absolute Ks restframe magnitude is of the order 
of -27.3 to -27.6. Such high luminosity objects should be ex- 
tremely rare; for this reason we subsequently only consider so- 
lutions with photometric redshifts below 4 and list the proper- 
ties estimated from the best fit models. 

A1835-#3: The best fit is obtained with the BC templates at 
Zflt = 1.12 for a burst of 0.5 Gyr with an extinction of Ay — 0.8. 
The corresponding stellar mass is M* ~ 5.1 x lO 9 /ji/M . Fits 
with the Maraston templates yield very similar parameters. A 
very similar mass (5.0 x 1O 9 //vM ) is also obtained using Eq. 
1 and L K /M = 3.2 adopted for SCUBA galaxies (Borys et al. 
2005). 

A1835-#10: The best fit is obtained with the BC templates 
at Zflt = 1 .23 for a burst of 0.5 Gyr with an extinction of Ay = 
1.8. The corresponding stellar mass is M* ~ 9.5 x 1O 9 /a/M , 
basically identical to the one derived following Borys et al.. 
Fits with the Maraston templates yield a somewhat younger 
age (0.25 Gyr), lower extinction (Ay = 0.4), and lower mass 
(6.4 x lO 9 /^M )atz fit = 1.14. 

To illustrate the expected IR to sub-mm SED of these ob- 
jects we show, in Fig.[T3] the best fits for #3 and # 10 obtained 
with the GRASIL templates. For reference the best fit redshifts 



obtained with the GRASIL templates are zm = 0.935 and 1.165 
for #3 and #10 respectively, and the SFR is ~ 2.4 and 4.3 M Q 
yr _1 , assuming fx = 1.4 for both sources. The global SED and 
IR to sub-mm fluxes of these objects are quite similar to those 
of #17. Remember, however, that the sources #3 and #10 are 
blended with neighbouring object in the Spitzer images, hence 
requiring high spatial resolution observations to potentially re- 
solve them with future instruments. 

5.6. The importance of IRAC/Spitzer photometry on 
the photometric redshifts 

For the objects with measurable IRAC photometry we have 
also examined the importance of this additional information on 
the SED fits and the resulting photometric redshifts. Overall 
this exercise, summarised in Fig. [T4l shows that the objects can 
be grouped in three "classes": 

1) Objects showing degenerate/ambiguous low- and high-z 
solutions even with IRAC photometry (cf. Fig. [3] with A1835- 
#2 shown in the top panel, also: A1835-#l to some extend). 
For A1835-#2, formally the best fit is found in both cases (i.e. 
with or without IRAC photometry) at z ~ 7.-7.5). However, the 
solution at lower z is clearly favoured from arguments on the 
absolute magnitude of this object (cf. Sect. [4j. 

2) Degenerate/ambiguous low- and high-z solutions whose 
degeneracy is lifted thanks to IRAC photometry (see A1835- 
#4, middle panel, and also AC114-#1). In this case the "cur- 
vature" measured between near-IR and the IRAC photometry 
allows one very clearly to locate the 1.6 yum stellar bump due 
to the minimum in the H opacity and hence to constrain the 
galaxy redshift. A1835-#l is intermediate between this case 
and the following one (3). 

3) Unconstrained photometric z from ground-based pho- 
tometry, which becomes well defined low-z solution with 
IRAC photometry (A1835-#17, bottom panel). This object, the 
faintest of our EROs in H and Ks, is a /-dropout. Therefore 
relying only on two ground-based near-IR photometric points 
results in a basically unconstrained photometric z, as shown in 
Fig [T4](left bottom panel, dashed line). The two IRAC detec- 
tions at 3.6 and 4.5 pm together with the upper limits at longer 
wavelength, are again attributed to the stellar 1.6 fim bump, 
making this object a very clear z ~ 0.81 galaxy. 

In conclusion, if sufficient near-IR ground-based and 
IRAC/Spitzer photometric datapoints can be secured and if 
they reveal a strong enough "curvature" such as to constrain 
the redshifted stellar 1 .6 //m bump, SED fitting of EROs can 
yield fairly well defined photometric redshifts. In other cases, 
ambiguities between low- and high-z may remain, and other ar- 
guments/data are needed to securely determine their redshift. 
However, given the resulting bright absolute magnitudes of 
these objects - if at high-z - and their large number we con- 
sider the "low-z" (0 ^ z ^ 3) solution much more likely for all 
of the objects studied here. 

Although not included in the SED fitting procedure done 
with Hyperz, a measurement or a stringent upper limit of the 
24 /jm flux with MIPS can also provide important constraints, 
helping e.g. to distinguish between SED fits with or without 
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strong redenning and therefore indirectly to rule out certain red- 
shifts (see e.g. the cases of A1835-#2, AC114-#1). 

6. SED fitting of related objects 

6.1. HUDF-J2 revisited 

Using the same modeling technique and spectral templates de- 
scribed above, we have examined the optical (zssolp) drop-out 
object HUDF-J2 discussed by Mobasher et al. (2005) and more 
recently by Dunlop et al. (2006). This object was also selected 
by Yan et al. (2004) on the basis of its red 3.6 pss\j z^soLP-band 
flux. Their full sample, named "IRAC" EROs (IEROs), will be 
discussed below (Sect. l6.2l i. 

6.1 .1 . Available photometry 

Here we have used two datasets published for this ob- 
ject, (i): JHK photometry from NICMOS and ISAAC, 
IRAC/Spitzer data (channel 1—4), and the zgsoLp non-detection 
from Mobasher et al. (2005). In addition the 24 pm flux mea- 
sured by Mobasher et al. (2005) is also used for comparison, 
though not included in the Hyperz spectral fitting, (ii): same 
data as (i), except for revised values of the optical data from 
Dunlop et al. (2006) yielding a B435 non-detection plus detec- 
tions in V606, / 77 5 and Z8 5 o LP . Where necessary, quoted non- 
detection limits have been transformed to 1 <x limits for a con- 
sistent treatment with the Hyperz photometric redshift code. 

6.1.2. Results 

Using the photometric dataset (i) from Mobasher et al. (2005) 
we find quite similar results as these authors, i.e. a degeneracy 
between solutions at z ~ 2.5 plus significant extinction and at 
z ~ 6-7.5 with zero/little extinction, with the formal best fits 
lying at high-z. This is illustrated by the;^ 2 maps shown in Fig. 
[TBI The main results for this object are summarised in Table [7] 
Concerning the high-z solutions we note the following. The 
best fit Zph ot depends somewhat on the spectral templates used; 
indeed using BCCWW or s04gyr templates we obtain zg t ~ 
7.39 Of 2 = 1.5) or zm ~ 6.5 Of 2 = 1.2). In both cases the 
best fit is obtained with zero extinction. With the same tem- 
plates groups the best fit "low-z" solutions are zg t =2.59 and 
A v = 1.8 0f 2 =2.3) for the BCCCWW templates, and z fit =2.42 
and Ay = 3.4 (y 2 =2.6) for the s04gyr templates. Test calcu- 
lations have shown, introducing a minimum error of ~ 0.1- 
0.15 mag already modifies considerably the^ 2 map leading to 
less well constrained solutions. Given uncertainties in the de- 
termination of measurement errors and uncertainties in match- 
ing photometry from different instruments (mainly NICMOS, 
ISAAC, and Spitzer) such errorbars may be more realistic than 
the small errors quoted by Mobasher et al. (2005). A compar- 
ison of the photometry of this object published by Yan et al. 
(2004, their object # 2) and the measurements of Mobasher et 
al. (2005), showing differences exceeding several cr in various 
bands, is also illustrative for this purpose. Given the relatively 
small differences in the x 2 , the resulting exceptionally bright 
magnitude of this object ( Af esl (Ks) ~ -27., cf. Table Ell, and 



using the same "spirit" as for our EROs, we would conserva- 
tively favour a "low-z" interpretation for HUDF-J2 on this ba- 
sis. Quantitatively the main weakness of the low-z fits is the 
slight excess predicted in the zssolp band with respect to the 
photometry from Mobasher et al., as shown in Fig. [16] Actually 
the flux predicted with the BCCWW template (black line) cor- 
responds to a 2<t detection; a somewhat larger flux is predicted 
with the GRASIL templates. 

Recently, Dunlop et al. (2006) have questioned the very 
deep ZgsoLp non-detection limit quoted by Mobasher et al. 

(2005) , and they have performed manual photometry of 
HUDF-J2 in the optical bands. Their faint detection in V606, 
i'775, and Z850LP alters the balance of the^ 2 behaviour between 
low-z and high-z favouring solutions at z&t - 2. 15 +0.3 (Dunlop 
et al. 2006). However, imposing a measurement in a prescribed 
aperture may not be appropriate. If we use the same constraint 
and the BC templates we obtain zgt = 2.48, Ay - 3.0, and a 
burst of 0.26 Gyr age, in good agreement with Dunlop et al. 

(2006) . In any case, the model fluxes shown in Fig. [16] in the 
Z850LP band and at shorter wavelengths, are nicely bracketed by 
the measurements of Dunlop et al. and Mobasher et al. 

As qualitatively the SED of HUDF-J2 resembles that of our 
EROs (in particular that of A1835-#2), and given that "low- 
z" fits indicate significant extinction, it is instructive to ex- 
plore also spectral templates of dusty objects. Indeed, using 
the GRASIL templates a reasonable fit (y 2 =4.8) to the observa- 
tions (i) are found with the M82 template and additional extinc- 
tion of Ay = 2.8 at Zfk ~ 1.8 (see red line in Fig.fTTb. Adopting 
the photometry from (ii), an excellent fit (y 2 =0.9) is found with 
the HR10 template and no additional extinction for zsx - 2.31 
(blue line in Fig. [P7I ). Interestingly, both fits reproduce quite 
well the 24 pm flux observed by MlPS/Spitzer (Mobasher et 
al. 2005), which was not included in our fit procedure. This 
lends further credit to the explanation of HUDF-J2 as a dust 
rich galaxy at z ~ 2.3. For comparison, Mobasher et al. do not 
fit the 24 pm emission, and invoke other components - possibly 
an obscured AGN - to explain this flux. And this measurement 
is not considered by Dunlop et al. (2006) in their analysis. As 
also shown by Fig. [17] the predicted far-IR to sub-mm flux of 
HUDF-J2 should be within reach of existing/future facilities. 
Detections in this spectral range should definitely be able to 
distinguish between dust-free high-z solutions and the "low-z" 
fits favoured here. 

From the two GRASIL fits shown in Fig. [17] with the M82 
and HR10 templates at ztit = 1-83 and 2.52 respectively we 
estimate a bolometric luminosity of L\, Q \ ~ (4.2 - 6.2) x 10 n L Q 
(in the LIRG range) corresponding to SFR ~ 72-107 M yr . 

6.2. IRAC selected EROs in the HUDF 

To compare our objects with the IRAC selected EROs in the 
HUDF from Yan et al. (2004) discussed in paper I, we have 
subjected them to the same quantitative analysis using Hyperz. 
The results from the SED fits are summarised in Table [8] Note 
that for simplicity we list fit results using Bruzual & Chariot 
templates only, although computations were also carried out 
using the other template groups. Furthermore, to account for 
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possible mismatches between the photometry from different in- 
struments (ACS, NICMOS, ISAAC, and IRAC) and to allow 
for other uncertainties in the error quantification, we adopt a 
minimum error of 0.1 mag in all filters. 

6.2.1. Properties of IRAC EROs 

As already noted by Yan et al., four objects have quite uncer- 
tain photometric redshifts. #1, 2, and 10 have quite degenerate 
X 2 over a large redshift range. #5 yields the worst fit, due to an 
apparent flux excess at 8 fim, which could be an indication for 
an AGN contribution. Therefore we follow Yan et al. and ex- 
clude these objects from the subsequent discussion of the aver- 
age properties of this sample. However, it is worth mentioning 
that # 2 corresponds to the HUDF-J2 object from Mobasher et 
al. (2005) discussed above. 

The main difference with the results of Yan et al. is that, 
once allowing for a reasonable minimum error of say 0. 1 mag 
and including a wide variety of star formation histories and 
varying extinction, we are able to obtain good fits to all ob- 
jects (except those already mentioned) with standard Bruzual 
& Chariot templates, i.e. there is no need for composite stellar 
populations as invoked by Yan et al. (2004). In consequence, 
our best fit ages and extinction differ systematically from their 
analysis, as will be discussed below. Otherwise quite similar 
properties are derived from our more complete quantitative 
analysis. 

More precisely we find best fit redshifts ranging from 0.6 
to 2.8 with a median (mean) of 1.6 (1.9). The redshifts dif- 
ference with Yan et al. has a median value of 6z/(l + z) ~ - 
0.07. In particular we obtain significantly lower redshifts for 
two objects, #13 and 14. However their photometric redshifts 
turn out to be quite uncertain and strongly dependent on the 
template set used 0. We therefore consider them as quite in- 
secure. Spectroscopic redshifts are available for 3 sources of 
the Yan et al. sample, #9, #13, and #17 (see Daddi et al. 2005 
and update by Maraston et al. 2006). Except for the uncertain 
object #13 just discussed, the agreement with our photometric 
redshifts is excellent (|<5z| < 0.05). 

If taken at face value our best fits yield the following av- 
erage properties for these IRAC selected EROs (cf. Table[8]): a 
median (average) extinction of Ay = 2.0 (2.2), a median (av- 
erage) age of 0.5 (1.3) Gyr, and stellar masses M+ ranging be- 
tween 10 9 and 5. x 10 11 Mo, with a median (average) mass of 
0.5 x 10 11 (1.3 x lO'^Mgjf 

Their absolute M%*\Ks) covers -20.5 to -24.3 in AB 
mags, with a median of -22.3, corresponding to a median of 
M r ^ ga (Ks) = -24.2. This is somewhat brighter than M* from 
2MASS (-23.53 cf. Kochanek et al. 2001) but similar to M* at 
Z ~ 0.6-2.5 from Bolzonella et al. (2002) and Kashikawa et al. 
(2003) who find M* e a (Ks) ~ -24.3 to -25.0. 

As already mentioned these average properties are in good 
agreement with those derived by Yan et al. (2004), except that 

3 Using the s04gyr (Maraston) templates we find best fits at zm = 
2.3 (2.1) and (0.7) 1.3 for #13 and 14 respectively. 

4 Note that our models assume a Miller-Scalo IMF from 0. 1 to 125 
Mo, whereas Yan et al. (2004) adopt a Chabrier (2003) IMF. 



we find good fits with templates computed assuming standard 
star formation histories, i.e. that we do not need to invoke other 
more arbitrary composite stellar populations. In consequence, 
the bulk of the fits we obtain correspond to younger stellar ages 
and hence to higher extinction than Yan et al.. As already men- 
tioned by these authors, observations at longer wavelengths (in- 
cluding MIPS imaging) could help to distinguish between these 
solutions; in fact, according to Yan (2006, private communica- 
tion) 9 of the 17 IERO have been detected at 24 yum with fluxes 
above 20 //Jy providing support to our interpretation. 

For the subsequent comparison we adopt the properties de- 
rived here as representative values for the IRAC selected EROs. 

6.2.2. Comparison with our ERO sample 

Empirically our optical drop-out EROs share many properties 
with the IRAC selected EROs (cf. paper I), in particular sim- 
ilar optical to near-IR and IRAC/Spitzer colours, and similar 
Z850LP, near-IR and IRAC/Spitzer magnitudes, after correcting 
for our median magnification factor corresponding to ~ 0.4 
mag. In contrast to our EROs, which so far are undetected in 
all bands shortward of ZgsoLP, the IRAC selected EROs of Yan 
et al. (2004) are detected in optical bands (BVIz), although at 
very faint magnitudes (e.g. V ~ 27 to 30.). As the depth of this 
data, taken from the Hubble Ultra Deep Field, is deeper than 
our optical photometry, this is not incompatible with the two 
samples being very similar. 

The similarity between our EROs and the IRAC selected 
EROs of Yan et al. (2004) is also evident from the comparison 
of the derived properties (cf. Tables [5] and [8]), showing similar 
redshifts, extinction, stellar ages, and stellar masses. 

7. Discussion 

A fact worth mentioning concerning the objects discussed here 
is that one of the sources, AC114-#1, fulfills the criteria com- 
monly used to select z ~ 6 galaxies as /-dropout. Indeed for 
this object (/814 _ Z85olp)ab > 2.16 and it is undetected at all 
wavelengths shorter than Zg50LP. Thus the commonly adopted 
i dropout criterion for ACS observations, (Ijjsyr— Z85olp)ab > 
1.5 (cf. Stanway et al. 2003, Dickinson et al. 2004) yield- 
ing (/775w-Z85olp)ab > 1-85 once transformed to the 7gi4 fil- 
ter from WFPC2 (Sirianni et al. 2005), is satisfied by this 
ERO. Furthermore, the flux expected in B and V filters from 
our SED fits is sufficiently low (Vf606w ^ 29.2 in AB mag) 
that it would remain undetected in the GOODS imaging with 
ACS (cf. Giavalisco et al. 2004). Therefore such a low red- 
shift object, with an intrinsic lensing corrected magnitude of 
Z850lpab ~ 26., could potentially contaminate /-dropout sam- 
ples. Apparently this is not the case for the sample of IRAC se- 
lected EROs from Yan et al. (2004), which if detected in zssolp 
are also detected at shorter wavelengths. 

Obviously, observations in several filters longward of Z85olp 
are able to eliminate such objects with very red SEDs out to 
the IR. For example, examining J - H Stanway et al. (2005) 
find that their iV^-dropout sample is consistent with unred- 
dened high-z starbursts. On the other hand, in a recent analysis 
combining ACS and Spitzer imaging of the GOODS field, Yan 
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Table 5. Derived/estimated properties for optical dropout ERO galaxies with near-IR and Spitzer detections. Listed are the 
object ID (col. 1), the photometric redshift estimate (col. 2), the extinction (col. 3), the type of the best fit template (col. 4), the 
distance modulus corresponding to z p hot (col. 5), the absolute Ks-band magnitude non-corrected for lensing (col. 6), the absolute 
rest-frame A's-band magnitude non-corrected for lensing (col. 7), the estimated stellar mass (from scaling the SED fit or from 
M^iKs) assuming L^/M = 3.2, col. 8), the estimated star formation rate non-corrected for lensing (col. 9), and the age of the 
stellar population (col. 10). To correct the above mentioned absolute quantities for gravitational magnification the appropriate 
magnification factors listed in Table |4]must be used. 



Object 


£phot 




A v 


template 


DM" 


M KSi (K s) - 2.5 logOu) 


MassXyu 


SFRXyU 


stellar age 








[mag] 




[mag] 


[mag] 


M e 


M Q yr 1 


[Gyr] 


A1835-#l 


-0.4 


-1.5 


? 


Fits uncertain - 


see text 










A1835-#2 


-2.8- 


-3 


2.4-3 


young burst 


47.0 


-27.7 


- 1.2 x 10 12 


-2100 


<0.36 


A1835-#4 


- 1.2 




0-1.6 


burst/elliptical 


44.60 


-21.6 


- 1.7 x 10 10 


-5 


0.7 to 4.5 


A1835-#17 


-0.7- 


-0.8 


- 3.8 


burst 


43.0 


-21.7 


- 1.3 x 10 10 ? 


-0.9 


? (see text) 


AC114-#1 


- 1.3- 


-1.6 


- 1.6-2.8 


burst 


44.84 


-26.4 


(1.3- 2.6) x 10 12 




- 0.9-4.5 Gyr 


AC114-#1 


- 1.0 




+3.8 


M51 


44.03 


-25.9 




-48 





" distance modulus computed for minimum redshift 



Table 6. Same as Table [5]for optical dropout ERO galaxies detected only in the near-IR (no Spitzer photometry available). No 
information is listed for A1835-#l 1 due to its highly uncertain photometric redshift. 



Object 


£phot A v 


template DM M lest (Ks) - 2.5 logOj) 


MassXyu 


SFRx/z 


stellar age 




[mag] 


[mag] [mag] 


M 


M yr" 1 


[Gyr] 



A1835-#3 -1.1 -0.6-0.8 burst 44.4 -22.2 - 5.1 x 10' 0.5 

A1835-#10 -1.2 -1.8 burst 44.68 -22.9 - 9.5 x 10 9 0.5 

A1835-#ll ? 



Table 7. Same as Table [5] for the HUDF-J2 galaxy from Mobasher et al. (2005). To the best of our knowledge no lensing 
correction has to be applied to this object (fj. — 1). 



Object 






A v 


template 


DM" 


M les, (Ks) 


Mass 


SFR 


stellar age 








[mag] 




[mag] 


[mag] 


M 


M Q yr 1 


[Gyr] 


HUDF-J2 


-2.4- 


-2.6 


- 1.8-3.4 


burst 


-46.47 


-24.1 to -24.5 


- 5. x 10" 




i 0.6-2. 


HUDF-J2 


-6.4- 


-7.4 


-0 


burst 


49.03 


-26.9 


(2. -3.)x 10 12 




S0.7 


HUDF-J2 


-2.5 






HR10 


46.58 


-23.9 




107 




HUDF-J2 


- 1.8 




- +2.8 


M82 


45.73 


-23.0 




72 





et al. (2006) find that - 15-21 % of the IRAC detected i 175W - 
dropouts have very high flux ratios between 3.6 /im and the 
Z850LP band. These amount to ~ 4 % of their total i dropout, i.e. 
Z ~ 6 sample. Our object, AC1 14-#1, belongs very likely to the 
same class of rare objects. 

8. Conclusions 

We have undertaken a detailed analysis of the stellar pop- 
ulations and the extinction of a sample of 8 lensed EROs 
in Abell 1835 and AC114 from Richard et al. (2006) and 
related objects from the Hubble Ultra Deep Field (HUDF). 
The analysis includes in particular one known SCUBA galaxy 
(SMMJ 14009+0252), the z ~ 6.5 post-starburst galaxy candi- 
date from Mobasher et al. (2005), and the IRAC selected EROs 
from Yan et al. (2004). 

Empirically, these objects share a very red overall SED, 
similar colours and magnitudes, and very faint or absent flux 
in optical bands. In particular most of our EROs, originally se- 
lected as very red (R - Ks > 5.6) optical drop-out objects by 



Richard et al., have been detected with ACS/HST in deep zssolp 
images, as discussed in Hempel et al. (2006, paper I). 

The ACS, VLT, and IRAC/Spitzer photometry has been 
taken from paper I. To determine photometric redshifts and to 
simultaneously constrain the stellar population and extinction 
properties of these objects, we have used an updated version 
of the Hyperz code from Bolzonella et al. (2000) including a 
large number of synthetic, semi-empirical and empirical spec- 
tral templates. 

The main results from the SED fitting are the following: 

- The SED analysis including near-IR plus IRAC photom- 
etry for 5 of our objects, shows in most cases degenerate 
solutions between "low-z" (z ~ 1-3) and high-z (z - 6- 
7) for their photometric redshifts. Although formally best 
fits are often found at high-z, their resulting bright abso- 
lute magnitudes, the number density of these objects, and 
in some cases Spitzer photometry or longer wavelength ob- 
servations, strongly suggest that all of these objects are at 
"low-z". 
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Table 8. Derived/estimated properties for IRAC selected EROs from Yan et al. (2004) from Hyperz fits with Bruzual & Chariot 
templates with different star formation histories, and computed assuming a minimum error of 0.1 mag in all bands. Note: all 
magnitudes are in the Vega system (M(Ks)y eg . d = M(Ks)& B - 1.871). 



Object 


ZYan 




template 


age 


A v 


DM 


M ms \Ks) 


Mass 










[Gyr] 


[mag] 


[mag] 


[mag] 


M 


1* 


3.6 


2.5 


burst 


0.7 


2.4 


46.60 


-25.85 


3.4e+ll 


2*a 


3.4 


2.1 


burst 


0.7 


2.8 


46.12 


-25.06 


1.7e+ll 


3 


2.9 


2.8 


t = 1 Gyr 


2.6 


2.0 


46.82 


-25.9 


3.5e+ll 


4 


2.7 


1.6 


t = 1 Gyr 


4.5 


2.0 


45.43 


-24.02 


7.7e+10 


5* 


2.8 


3.0 


t = 1 Gyr 


2.3 


2.2 


47.04 


-26.22 


2.1e+ll 


6 


2.3 


1.1 


burst 


0.01 


3.8 


44.40 


-22.8 


1.2e+09 


7 


2.7 


2.6 


t = 5 Gyr 


2.6 


1.8 


46.69 


-25.66 


5.2e+ll 


8 


2.9 


2.7 


t = 1 Gyr 


2.3 


1.4 


46.71 


-26.13 


3.8e+ll 


9 


2.8 


2.7 


burst 


0.4 


1.0 


46.76 


-25.19 


1.4e+ll 


10* 


2.1 


0.9 


burst 


0.2 


3.8 


43.75 


-20.58 


8.3e+08 


11 


2.4 


1.6 


burst 


0.006 


3.8 


45.36 


-23.4 


6.2e+09 


12 


1.9 


2.5 


burst 


0.4 


0.6 


46.55 


-23.54 


2.9e+10 


13* 


1.9 


0.7 


burst 


0.2 


3.8 


43.00 


-22.38 


9.5e+08 


14'' 


1.6 


0.6 


t = 5 Gyr 


2.6 


3.8 


42.61 


-22.79 


2.1e+10 


15 


2.7 


2.7 


t = 1 Gyr 


1.7 


1.0 


46.75 


-24.16 


5.2e+10 


16 


2.4 


1.6 


burst 


0.01 


2.4 


45.44 


-24.33 


5.0e+09 


17 


1.6 


1.7 


burst 


0.5 


0.8 


45.50 


-24.38 


5.1e+10 



* Uncertain photometric redshifts 

" same as object HUDF-J2 from Mobasher et al. (2005) 

b Zphoi redshift likely underestimated (see text) 



- The majority of our lensed objects are best understood by 
relatively young (5= 0.5-0.7 Gyr) and dusty starbursts. 

- For 3 of our objects we find indications for strong ex- 
tinction, with Ay ~ 2.4-4. Among them, the galaxy 
SMMJ 14009+0252 is a known sub-mm emitter and at least 
one of them, the ERO AC114-# 1, is predicted to be in 
the LIRG category and expected to be detectable with cur- 
rent sub-mm instruments. For the remaining objects, among 
which 3 show moderate to strong extinction, we present 
predictions of their IR to sub-mm SEDs for future obser- 
vations with APEX, Herschel, and ALMA. 

- The stellar masses estimated for our objects span a large 
range from ~ 5. x 10 9 /// to 1. x lO 12 /yuM , where fi is the 
magnification factor derived from the gravitational lensing 
model. Typically one has y. ~ 1.2, with a maximum mag- 
nification of ^ 2. Where appropriate, star formation rates 
estimated from the bolometric luminosity determined from 
spectral template fitting are SFR ~ (1 — 36)/yu M Q yr _1 . 
For SMMJ 14009+0252, the most extreme case, we esti- 
mate Lboi ~ 6.xlO n L and SFR ~ 1000 M yr 1 forz fit ~ 3 
and /i = 2. 

- Taking uncertainties and revisions of the photometry of 
HUDF-J2 as well as accounting for a variety of spectral 
templates, we suggest that this object, originally identified 
as a z ~ 6.5 massive post-starburst, is more likely a dusty 
starburst at z ~ 2.3-2.6 in agreement with the reanalysis 
by Dunlop et al. (2006). We show that this explanation also 
naturally explains the observed 24 fim emission from this 
object and we predict its IR to sub-mm SED. 

- Using the same methods we have analysed the sample of 
IRAC selected EROs from Yan et al. (2004). Both empiri- 
cally and from our SED fits we find that these objects show 



very similar properties to our lensed EROs. In contrast to 
Yan et al. we do not find strong indications for composite 
stellar populations in these objects. 
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Fig. 1. Observed SED of all EROs in the optical and near-IR up Fig. 3. A1835-#2: Comparison of best fit high-z solution 

to 10 /im. The two most extreme objects, the sub-mm galaxy (black, Zfi t =7.46) with model fits at z ~ 3. Blue: Bruzual & 

A1835-#2 and AC114-#1, are indicated in red and yellow re- Chariot model (SF history as ellipticals at age of 0.36 Gyr + 

spectively. A1835-#l 1, whose properties including the photo- 2.4 mag A v extinction). Red: GRASIL template of NGC 6090 



metric redshift remain highly uncertain is shown in blue. 



lllER0aug06_bccww_matrix. 1 093. fits 




Fig. 2. x ma P as a function of redshift and extinction 
(x 2 (z,Av)) for A1835-#2, displaying in dark the most proba- 
ble regions on a logarithmic scale. Solid lines enclose the 1 
to 3 cr contours (confidence levels of 68, 90 and 99 % respec- 
tively). The (z,Ay) projection plane presented in this figure cor- 
responds to the best x 2 found through the model-age param- 
eter space for templates from the BCCWW group. Note the 
degeneracy of the photometric redshift solutions in this plane. 
Formally the best fit is found at high redshift (zg t > 7). However, 
for various reasons, including the 24 fim and sub-mm SED and 
the exceptional luminosity of this object if at high z, the most 
likely redshift of this galaxy is ~ 3. See discussion in text. 



with A v =1.4. See discussion in text. 



A1835#2: z = 7.46 (black), 2.95 (blue), 2.78 (red), 2.88 (green] 
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Fig. 4. Observed SED of the ERO/sub-mm galaxy A1835-#2 
including VLT, Spitzer (IRAC, MIPS), SCUBA observations. 
The model fits are the same as shown in Fig. [3] Note that the 
two fits with GRASIL models at z ~ 2.8-2.9 reproduce well 
the MIPS 24 yum flux and bracket the observed sub-mm points. 
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Fig.5. Same as Fig.|2]for A1835-#l. 

A1B35 #1 
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Fig. 8. A1835-#4: Comparison between fits with BCCWW 
galaxy template at z — 1.20, and fits with the M82 template 
from GRASIL at z = 1.24 (and no additional extinction). 

QllER0aug06_ulnov06_bccww_matrix,31 1 .fits 




1 2 
log(wavelength) [^m] 



Fig. 6. A1835-#l: Comparison between fits with Bruzual and 2 4 

Chariot models at z = 0.50 (black) and 1.35 (blue), and fits with 

dusty starburst models from GRASIL templates at z = 1.75 Fig.9. Same as Fig.|2]for A1835-#17. 
(M82, red) and 0.4 (green, M82 plus additional extinction of 
A v = 3.8. 

allEROaug06_ulnov06_bccww_rnatrix. 504. fits 




Fig. 7. Same as Fig.|2]for A1835-#4. 
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non-detected lensed EROs and related objects 
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Fig. 10. A1835-#17: Comparison between fits with Maraston 
(2005) templates at z = 0.81 and A v = 1., and fits with dusty 
starburst models from GRASIL templates at z = 0.60 and 0.79 
plus additional extinction of A v ~ 3.6. Note that the Maraston 
template does not include dust emission, which would be ex- 
pected for such high an extinction. 
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Fig. 11. Same as Fig.|2]for AC1 14-#1. 



Fig. 12. AC1 14-#1: Best fit SEDs with Bruzual & Chariot tem- 
plates (black line: burst of 4.5 Gyr age and Ay = 2.4 at z = 1 .3), 
S04gyr templates (red: 1.0 Gyr, A v = 2.8, z = 1.6), and with 
GRASIL templates (blue: M52 template + A v = 3.8, z = 1.0). 
Note that by construction only the GRASIL templates include 
dust emission. 
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Fig. 13. Model fits for A1835 #3 (black) and #10 (blue) using 
GRASIL templates showing predictions for the Herschel/sub- 
mm/ALMA spectral domain. 
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Fig. 14. Comparison of photometric redshifts derived with or 
without IRAC photometry for selected objects. Left panels: x 2 
of the best fitting template versus redshift obtained from SED 
fits including all the photometry (black solid lines), or obtained 
from SED fits excluding the IRAC photometric measurements 
(red dashed lines). Right panels: Best fit SED including all 
photometry (black solid curve) and excluding IRAC photom- 
etry (red solid). The observations are show by the blue sym- 
bols. From top to bottom the objects are: A1835-#2, #4, and 
#17. These objects illustrate three typical behaviours: 1) de- 
generate/ambiguous low- and high-z solutions even with IRAC 
photometry (#2, top panels), 2) degenerate/ambiguouslow- and 
high-z solutions whose degeneracy is lifted thanks to IRAC 
photometry (#4, middle panel), and 3) unconstrained photo- 
metric z from ground-based photometry, which becomes well 
defined low-z solution with IRAC photometry (#17, bottom 
panel). 



Fig. 16. z ~ 1.8-2.6 SED fits to the observations of HUDF-J2 
from Mobasher et al. (2005). Black line: fit with BCCWW tem- 
plates at z fit =2.59. Red: GRASIL M82 template and additional 
A v = 2.8 at Zfl, ~ 1.8. Blue: GRASIL HR10 template without 
additional extinction at ztit ~ 2.3. 
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redshift 

Fig. 15. Same as Fig. [2] for HUDF-J2. Formally the best fit is 
found at high redshift (zfit~ 6.4-7.4). However, for various rea- 
sons, including the 24 jjm flux and the exceptional luminosity 
of this object if at high z, we consider it more likely that this 
galaxy is at z ~ 2.4-2.6. See discussion in text. 
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Fig. 17. SED fits to the observations of HUDF-J2. The near- 
IR and IRAC/Spitzer photometry is taken from Mobasher et 
al. (2005); in the optical we adopt either ((i) the zssolp non- 
detection from Mobasher et al. (2005) or ((H) the revisited 
photometry from Dunlop et al. (2006) yielding a B435 non- 
detection plus detections in V606, '775 and zssolp- Fitting ((i) 
with BCCWW templates yields a best fit at high-z (z ~ 7.4, 
black line). Fitting ((i) with GRASIL templates, we obtain 
a best fit at z ~ 1.8 with the M82 template and additional 
A v = 2.8 (red line) or at z ~ 2.3 with the HR10 template (and 
no additional extinction; shown in blue). Adopting (ii) the best 
fit with GRASIL templates is at z ~ 2.3 with the HR10 tem- 
plate (and no additional extinction; shown in green). Note that 
all "low-z" solutions using GRASIL templates reproduce natu- 
rally the observed 24 jjm flux. 



